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Effects caused by the event-by-event ffuctuation of the initial conditions in hydrody- 
namical description of high-energy heavy-ion collisions are investigated. Non-negligible 
effects appear for several observable quantities, even for a fixed impact parameter b . They 
are sensitive to the equation of state, being the dispersions of the observable quantities in 
general smaller when the QGP phase appears at the beginning of hydrodynamic evolution 
than when the fluid remains hadron gas during whole the evolution. 

1. INTRODUCTION 

In usual hydrodynamic description of high-energy heavy-ion collisions, one customarily 
assumes some highly symmetric and smooth initial conditions, which correspond to mean 
distributions of velocity, temperature, energy density, etc., averaged over several events. 
However, our systems are not large enough, so large fluctuations are expected. What are 
the effects of the event-by-event fluctuation of the initial conditions? Are they sizable? 
Do they depend on the equation of state? Which are the most sensitive variables? These 
are some questions which arise regarding such an initial-state fluctuation, and we try to 
shed some light on these matters in the present study||I|. 



2. METHOD OF STUDY 

In order to study the problem stated above, first we generate events by using the NeXus 
event generator 0, from which initial conditions are computed at the time r = 1 fm. Then, 
the hydrodynamic equations are solved, starting from these initial conditions, assuming 
some equation of state (EoS). To see the EoS dependence of the effects we are treating, 
we consider two different EoS's0: 

1. Resonance Gas (RG): c^^ = 0.2; 



*Work supported in part by FAPESP (contract nos. 2000/04422-7 and 98/00317-2), FAPERJ (contract 
no.E-26/150.942/99), PRONEX (contract no. 41.96.0886.00) and CNPq-Brasil. 



2 



0.2 , e< 0.28 GeV/W, 

2. QGP+RG: Cs'^ = { 0.056/e, mixed phase, 

1/3 - 45/3e, e > 1.45 GeV/W. 

The resolution of the hydro dynamic equations deserves some special care, since our 
initial conditions do not have any symmetry nor they are smooth. We adopt the so-called 
smoothed-particle hydrodynamic (SPH) approach^], first used in astrophysics and which 
we have previously adapted for heavy-ion collisions|0, a method flexible enough, giving 
a desired precision. The main characteristic of SPH is the parametrization of the flow in 
terms of discrete Lagrangian coordinates attached to small volumes (called "particles") 
with some conserved quantity. In the present work, besides the energy and momentum, we 
took the entropy as our conserved quantity. Then, its density (in the space-fixed frame) 
is parametrized as 

N 

s*{^,t) = J2u, wi^-Mty,h) , (1) 

i 

where 

iy(x — Xj(t); /;,) is a normalized kernel; 

Xj(t) is the i-th particle position, so the velocity is Vj = d'^.-Jdi ; 
h is the smoothing scale parameter; 

and we have 

N 

S = Jd^^s*{^,t) = Y.Ui. (2) 

i 

The equations of motion are then written as the coupled equations 



d / Pi + Si 



— V. 



dt 



)+Y.^J + V.l^(x, - X,; /i) = . (3) 



Following this procedure, we computed some observable quantities, event-by-event, for 
y/s =130y4GeV Au + Au collisions. The results are presented in the next Section. 

3. RESULTS 

3.1. Elliptic flow coefficient t>2 

Having solved the coupled equations (3), we have computed the particle spectra at 
T = and from which the elliptic flow coefficient V2 on an event-by-event basis. In 
Figure 1, we show its distributions for a fixed impact parameter b, for the two EoS 
considered. As expected, V2 exhibits a large fiuctuation, which depends on the EoS. 
One should take care in looking at this Figure that our b is the true impact parameter 
(not determined in the way experimentalists do), so for instance in the RG case, there 
are some events with negative V2 , which experimentally would not appear. As for the 
average values < f 2 > ? it is almost independent of the EoS. This is shown in Figure 2, 
where < f 2 > ± Sv2 is plotted as function of the centrality and compared with data[||. It 
is seen that < f 2 > reproduces well the experimental trend, whereas the dispersions 6v2 
are much wider than the experimental errors. As for the EoS dependence, (5f 2 is smaller 
when QGP is produced. 
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Figure 1. Distribution of elliptic-flow coefficients f 2 at 6 = 10 fm for two EoS. 
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Figure 2. EoS dependence of < f 2 > ±5f2 Figure 3. Pion m-r spectra for 5 events 

as function of the centrality, compared with [b = 7.0 fm). 
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3.2. rriT distributions 

In Figure 3, we show the ttit distributions for 5 events. As expected, mj- distributions 
are in general steeper when QGP is produced. As for fluctuations, the resultant fluctuation 
in rriT spectrum (or in the slope parameter ST) is very small. 

3.3. Multiplicity fluctuation in the mid-rapidity region 

Table 1 summarizes the results of our study on multiplicity fluctuation in the mid- 
rapidity region. It is seen that as 6 — > 0, < n^r > becomes much larger with the QGP 
EoS; ii) shows the same tendency in this limit; in) As for the ratio (5n^/ < >, it 
is not sensitive to the EoS. 
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Table 1 

EoS dependence of the multiplicity fluctuation in two different rapidity intervals —Ay < 
y < +Ay, as function of the impact parameter b. 



# of 

b [fm] EoS 

' events 


Ay=1.875 


Ay=3.00 


<n> 6n Sn/<n> 


<n> Sn 6n/<n> 


RG 44 
QGP 38 


1029.7 46.2 0.045 
1553.0 80.9 0.052 


1623.2 68.7 0.042 
2544.6 129.0 0.051 


RG 55 

7 

QGP 58 


613.3 49.5 0.081 
926.1 81.1 0.087 


977.4 71.6 0.073 
1530.5 123.7 0.081 


RG 166 
QGP 180 


312.8 43.0 0.137 
437.5 66.5 0.151 


506.1 65.8 0.130 
740.7 103.9 0.140 


RG 79 

1 2 

QGP 100 


162.8 35.6 0.219 
220.1 52.8 0.240 


268.9 56.2 0.209 
379.8 85.2 0.224 



4. CONCLUSIONS AND OUTLOOK 

The present study shows that the effects of the event-by-event fluctuation of the initial 
conditions in hydrodynamics are sizable and should be considered in data analyses. They 
do depend on the equation of state. Among the quantities examined here, 6v2 is the most 
sensitive to the equation of state. 

In the present work, many important factors have not been considered: baryon-number 
conservation, strangeness production, resonance decays, continuous emission effects, spec- 
tators, etc., which should indeed taken into account in order to get more precise results. 
Especially, use of the same procedure for the determination of the centrality as used by 
experimentalists, as in0, will make the results more directly comparable with data. In 
any event, we believe that the effects we studied will be present and will be sizable, even 
with these improvements. 
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